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1. Introduction 

A recent paper |I| demonstrated the existence of a family of closed-form solutions to 
the one-dimensional Schrodinger equation for the potential 

Vi{x) = -^4t^ (1) 

which is widely used 0, ^, |^ in the modeling of atomic response to strong time- 

dependent radiation fields. The approach used in Ref. was somewhat analogous to 
the Sturmian method p for generating a complete basis set of Coulomb wavefunctions, 
in which the energy E is held fixed and the nuclear charge Z is determined as an 
eigenvalue. Ref. ^ determined "eigenvalues" a„ of the cutoff parameter a, for which 
the eigenenergies En of the Schrodinger equation for the potential Eq. ([l|) take the 
values 

= 

these are the energies of the states of the three- dimensional hydrogen atom, in the usual 
system of atomic units (used throughout the present paper) in which the numerical 
value of unity is assigned to the mass m and charge e of the electron and to the reduced 
Planck's constant h. The associated eigenfunctions "^ni^) were found to have a closed- 
form expression: 

vi>„(x) =x^e-'^>^^(x), (3) 



where: z/ = or 1 for cases of even and odd parity respectively; x = + a„^; /„ is a 
polynomial of degree n > 1 — u; k, = 1/n; and root of an n -degree polynomial. 

In the present work a similar approach is applied to the three-dimensional smoothed 
Coulomb potential 

= - . ,f (4) 



This potential arises in the Kramers-Hennenberger transformation of the equations of 



motion of a hydrogen atom in a radiation field |]I0|, |Il|] , and it is qualitatively similar to 
some pseudopotentials used in density-functional calculations of electronic structure (e.g. 
Ref. |12[)- It is shown here that, subject to an assumption similar to one made in Ref. 



the Schrodinger equation for this potential admits an infinite number of closed- form 
solutions for each value of the angular momentum /. The lowest-energy solutions for 
each / are quite simple, so knowledge of these exact results may be useful for calibrating 
numerical methods that must be used to solve general equations of this type. To the best 
of the author's knowledge, these are the first exact results for Schrodinger eigenvalues 
and eigenfunctions for potenials of the class described by Eq. (H). 
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2. Method of solution 

The extension of the method of ||T[| is straightforward. We proceed from the Schrodinger 
equation 

-W^ir) - vl/(r) = E^ir), (5) 

and invoke the usual spherical coordinate factorization \l'(r) = ilJi{r)Yim{0, 0). With the 
definitions a = Za, p = ^ [ZrY + a^, e = E /Z"^ = — and 

Hr)=r'Up). (6) 

the equation 

is obtained. We now postulate a solution of the form 

n 

0/(p) = e-''''Ec.p\ (8) 

i=0 

where n is an integer to be determined. For economy of expression the dependence of k 
and Cj upon n and / is not made explicit in the notation, but it is held to be implicit. 

We now show that Eq. does not provide a general solution, but that it is 
applicable for discrete values of a. To see this, we substitute eq. (|]) in the l.h.s. of 
Eq. (^), and require the net coefficient of each power of p to vanish identically. The 
coefficient of fP that emerges from this operation is 



-a\3 + 4)(j + 2)c,+4 + «'/€(2j + 5)c,+3 + 

((j + 2)(j + 2/ + 3) - o'k') + 2 (1 - k(j + / + 2)) c,+i = (9) 

for i = —3, —2, — 1, n — 1. Thus we must solve the n + 3 simultaneous equations (9) 
in the n + 2 variables {ci, C2, c„, k, a}. 

We proceed by treating Eq. (9) as a four-term recurrence relation, and solve for 
Cj sequentially downwards from j = n — 1. Because cj = for j > n, the solution for 
j = n — 1 gives 

K = ^ . (10) 

Thus K is uniquely determined by n and I. Eq. (|T^) describes exactly the spectrum of 
the three-dimensional Coulomb potential, i.e. the limiting case of Eq. with a = 0; in 
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that case n coresponds to the number of nodes in the radial eigenfunction. We shall see 
subsequently that, for a 7^ 0, we must have n > 1, and the maximum number of radial 
nodes in the wavefunction described by Eqs. (||) and (|I0|) is n — 1. The potential of Eq. 



(^) has a long-range Coulomb tail and a non-Coulombic component at small r, so its 
Schrodinger spectrum is naturally described in the language of quantum defect theory 
I3| . In that terminology, Eq. (|I0|) describes a state with an integer value of the quantum 



defect /i, which is necessarily negative for a 7^ 0. The following development will indicate 
that all eigenfunctions with integral quantum defects obtained with potentials of the 
class Eq. (|) are described by Eq. (|). 

Since Eq. is homogeneous, we can set c„ = 1 without loss of generality. The 
values of Cj for i < n are then determined in terms of c„ by downward recursion using 
Eq. (9). For j = n — 2, we get 

c„_i = I (n{n + I + l){n + 2/ + 1) ^^-^ ) c„ = Pi(a^)c„, (11) 

2 \ n + I + 1J 

where pi{x) designates a first-degree polynomial in x. From this equation it is apparent 
that there will be a solution for n = only if a = 0, which is the familiar Coulombic 
case. 

Inspection of the structure of Eq. (9) shows that by continuing the recursion process 
downward in j we get 

with Pm(«^) being a polynomial of degree m in ct^. So for j = and j = — 1 we find 

ci = p„_i(a^)c„ and Cq = p„(a^)c„ (13) 

respectively. The last two cases to be considered are j = —2 and j = —3. These are 
found to give the same equation: 

^1 = = , ? , 1 ■ (14) 
n + I + 1 



Thus, from Eqs. ([131) and (|T4[)) we see that Eq. (^ will provide a valid solution if 

gn(a') = Pn(a') - (n + / + l)p„-i(a') = (15) 

or, in other words, if is a root of the n*'* degree polynomial g„. 

The applicability of Eq. (||) thus depends upon some roots of Eq. (|1^) being positive 
real numbers. The investigation reported here has not uncovered a general proof that 
Eq. (|15|) has any such roots, but calculations carried out for / < 10^ and n < 20 suggest 
that all its roots are positive real numbers and are nondegenerate. Let us adopt this as 
a hypothesis. If it is true, then the following statements hold: 
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i. For each l,n there are n + 1 values of o? for which Eq. ij^) has solutions of the 
form Eqs. (0, |^. This includes the previously-known (Coulombic) value = 0, plus 
the n roots of the polynomial g„. 

ii. The only potentials of the class eq. which have bound states with Coulombic 
energies are just those with values of that are solutions to Eq. ([73|j. This is because 
Eq. describes a monotonic function of a^. Thus, for a given /, its associated discrete 
Schrodinger eigenvalues will increase uniformly towards zero as a increases. For a given 
Z, Z, and Coulombic energy, 

Em = -\ . ,2 , (16) 

^ (n + / + 1) 

there will be some maximum value of o? for which E^i occurs as an eigenvalue, i.e., that 
for which it is the lowest eigenvalue. As cP' is decreased from this maximum, E^i will 
next occur in the spectrum when it is the second lowest eigenvalue, then as the third 
lowest, etc., until finally at o? = 0, when it is the {n + l)-th lowest eigenvalue. Thus 
there are indeed only n + 1 values of for which Eq. ([TBI) occurs in the spectrum, 
which is consistent with the stated hypothesis, and if the hypothesis is true, these values 
must thus coincide with the roots of Eq. (|T5|) 

The argument ii. is illustrative of the actual results of computations of solutions 
of Eqs. (9) and ([l5|), as will be described below. The largest value of for a given 



(n, /) is associated with a nodeless eigenfunction; the next largest value of a^, with 
an eigenfunction with one node; and so on to the smallest nonzero value of a^, which 
corresponds to an eigenfunction with n — 1 nodes. 

3. Results 

The eigenvalues of can be easily be found for a given (n, /) by numerical solution of 
the polynomial equation (|T5|) . They have relatively simple closed forms for n = 1 and 2, 
which are presented here. Numerical tables are given below for 1 < / < 3 and n < 10. 

3.1. n = l 

For n = 1 we obtain 

a2 = 2(/ + 2)^ Co = 1 + 2, ci = l, (17) 
so that the (nonnormalized) solution of Eq. (|^) is 

^^(r) = r'e-^/('+2)(/ + 2 + p), (18) 



with p = J{Zrf + 2{l + 2)3, and = Z-'^a^ = Z-^2{1 + 2^ . 
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3.2. n = 2 



For n = 2 there are two solutions for a^: 

(/ + 3)l 



a 



3± 



// + 15 



/ + 3 



(19) 



With the choice of C2 = 1, we can write the expressions for the two sets of coefficients 
Ci in the common form: 



Co 



- (/ + 3)2(2/ + 3), 



Cl 



Co 

/ + 3' 



(20) 



with the value of to be chosen as appropriate. 



3.3. n>2 

Although closed-form expressions can be obtained for and q for n = 3 and 4, they 
have the typical cumbersome form of roots of cubic and quartic equations, and it does not 
seem particularly useful to record them here in full. However, a simplifying relationship 
is worth noting. For these values of n, the values of can be written as 



a 



where /? is a root of the polynomials 



4^/5 



/ + 4 



/ + 14 
24-—— = 



(21) 



(22) 



for n = 3 and 4, respectively. In the limit of large /, the solutions of Eqs. (|2^) and (p^) 
tend respectively to /5 = 0, ±2 and P = ±1, ±3. Thus from Eqs. ( p!7| - |23|) we see that for 
large /, the smallest value of tends to = 2{n + / + 1)^ for n = 1 through 4. This 
motivates the choice of Eq. (|2TD as a general representation for the values of a^, and it 
has been used to record those values in the tables given below. It has been found that, 
to a high degree of numerical accuracy, the computed values of j3 for a given (n, I) sum 
to zero, so that substituting /5 = in eq. ^ apparently locates = (n + / + lY{n + 1) 
as the average of the values of a^. No fundamental explanation of this apparent fact is 
advanced here. 

Tables |I], 0, | give the values of /? for n < 10 for / = 1 — 3. A similar set of values 
for / = can be found in Table 1 of Ref. so they are not repeated here []T1 . 
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3.4- Systematic behavior of wavefunctions 

Numerical calculations indicate that the wavefunctions described by Eq. exhibit the 
qualitative behavior discussed at the end of Sec. For a given (n, /), denote by 
the /c*'^ smallest value of a obtained in solving Eq. ([15|) , with = 1,2, ...,n. Numerical 
experiments indicate that the wavefunction corresponding to ak has rij. = n — k radial 
nodes, a pattern that was observed in the one-dimensional cases treated in ref. [|l|. An 
example of this behavior is illustrated in Fig. for the case n = 10, / = 1. 

Fig. I depicts the values of afc for < Z < 4 as a function of ri,. + / + 2 = n*, 
the effective principal quantum number of atomic spectroscopy, which is related to the 
energy via Eq. (|ToD. For a given {l,k), as labelled in Fig. ||, a discernable sequence 
of values of ak is observed; these sequences are seen to approach definite limits as n* 
increases. This is related to the well-known phenomenon in atomic spectroscopy in which 
quantum defects tend to constant values high in Rydberg series. The slow variations of 
high-TT, quantum defects are due to the presence of a fixed, short-range, non-Coulombic 
part of the potential experienced by a Rydberg electron. Correspondingly, the fixed 
value of quantum defect obtained in the present method is associated with slow variation 
of ttfc at large n. 

For large and large k, on the other hand, the eigenfunctions approach those of the 
three-dimensional harmonic oscillator. This can be seen by expanding Z/ v^rM-~a^ in 
powers of r/a for large a; retaining the lowest two terms gives the Schrodinger equation 
for the harmonic oscillator. If first-order perturbation theory is used to include the 
effects of the term in this expansion, we obtain the approximate spectrum, 

^ + [2n, + / + 3/2] 

-A [6^r K + / + 3/2) + (/ + 5/2) (/ + 3/2)] , (24) 

as a — s> oo. Fig. ^ is a correlation diagram that displays the connection between this 
limit and the hydrogenic limit a = 0, if one keeps the number of radial nodes, rir, 
fixed as a varies. Two familiar cases of /-degeneracy are apparent in this figure: + I 
= constant for the hydrogen atom, and 2np + I = constant for the three-dimensional 
harmonic oscillator. 

3. 5. The spectrum for large values of I and a 

In sec. |3.3| it was mentioned that as / — oo, we find a„ 2{n + / + 1)^ for 1 < n < 4. 
Thus in this limit we recover a case of near-/-degeneracy similar to that encountered in 
hydrogen: there are values of a that support degenerate eigenfunctions with different 
values of /, described as a class by the equation + / = constant. The approach to this 
limit is relatively slow, apparently like as suggested by Eq. (|T^: e.g. for / = 10, 000 
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and n < 10, the actual value of a changes by about a part in 10^ for a unit change in 
at constant + 1- Numerical experiments suggest that this approximate degeneracy 
is a general phenomenon at large /. 

There is a simple effect of this kind of degeneracy for all potentials that have a long- 
range Coulomb tail and some non-Coulombic behavior localized at small r. In such 
systems, the centrifugal barrier presented to a high-/ wavefunction - commonly called 
a "nonpenetrating orbital" - will prevent it from sampling the non-Coulombic region. 
Thus the Coulombic /-degeneracy is largely undisturbed for large /. The quantum defects 
of nonpenetrating orbitals tend to zero as / ^ oo, a phenomenon that is universal 
in actual atomic systems, where, apart from isolated instances of series perturbation, 
observed quantum defects are hardly greater than 0.01 for / > 5 . 

However, the effect encountered in the present system is quite different. The non- 
Coulombic behavior of the potential extends to very large r, so the large-/ eigenf unctions 
are substantially modified from their Coulombic forms: their quantum defects are 
negative integers. The appearance of this novel /-degeneracy presumably derives from 
the existence of a constant of motion for the Schrodinger equation given by Eq. (^) that 
emerges in the large-/ limit, but it has not been identified in the present work. 

4. Conclusions 

A simple algebraic method has been presented to generate an infinite number of 
parameters a for which closed-form solutions may be found to the Schrodinger equation 
for the class of smoothed Coulomb potentials described by Eq. (^). The procedure 
bears some superficial resemblances to the Sturmian approach. However, it is not 
based on a system of orthogonal polynomials, and because the functions it generates 
are obviously incomplete, it probably cannot be simply related to known orthogonal 
systems. These results should be useful for testing, to arbitrary numerical accuracy, 
methods that integrate the Schrodinger equation for Coulomb-like systems, such as are 
encountered in electronic structure and collision problems. The approach also points to 
the possibility of previously unknown integrals of motion in these systems. 
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Figure captions 



Figure 1. Solutions for n = 10, I — 1: wavefunctions ipi{p) (solid line) and 
scaled potentials Z^V{r) — — [p^ + a^] ^ (dashed line) vs. p (horizontal axis). The 
numerical value of a is displayed for each case. The horizontal range < p < 350 is 
the same for all figures, and the negative portion of the vertical axis that is displayed 
covers the range [-0.11,0] in all cases. The wavefunctions ipi{p) are not normalized, and 
have each been scaled to fit the frame; they have also been chosen to be positive near 
p = 0, a, convention that differs trivially from that used in the constructive procedure 
presented in sec. ^, where instead the long-range tail is taken to be positive. 



Figure 2. Plot of a vs. n* = + ^ + 2 for < ? < 4, 2 < n + ? + 1 < 16, and a < 100. 
Key: plus: / = 0; x: / = 1; star: I = 2; square: / — 3; triangle: I = 4. The first few 
series {l,k) are exphcitly labelled. 



Figure 3. Correlation diagram for the spectrum of eq. (||) as a function of a (increasing 
schematically to the right), refiecting the conservation of n^. The vertical position of a 
level is proportional to its effective principal quantum number n*, which determines the 
energy via E = -Z^/(2n*^) or = -Z/a+ {n* + 3/2)y/Z/a^ for the hydrogenic and 
oscillator limits, respectively; the horizontal position corresponds to I as indicated. The 
values of n* for the two limits arc displaced and on different scales to ease visualization 
of the reordering of levels. 
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Table 1. Values of (3 as defined in Eq. (|i|), for Z = 1, 3 < n < 10. The n values of (3 
for each n are given in ascending order in the columns as labeled. 



n 



6 



-3.436527918374 
-0.954320351864 
4.390848270238 



-4.654288835724 
-3.021850153582 
0.696058817445 
6.980080171861 



-5.773876395778 
-4.650636370300 
-1.972562395354 
2.712910374971 
9.684164786461 



-6.844547165256 
-6.042758624724 
-4.058664349356 
-0.480918834922 
4.966113875958 
12.460775098300 



n = 7 



9 



10 



-7.888799456947 
-7.298711510457 
-5.793704911519 
-3.007931174911 
1.322982663950 
7.380032760853 
15.286131629031 



-8.917837293622 
-8.472195004982 
-7.307257906333 
-5.101818236107 
-1.607950021737 
3.353214980339 
9.908007968597 
18.145835513845 



-9.937641169524 
-9.593574030219 
-8.675819112578 
-6.904370845863 
-4.050370418227 
0.059395253271 
5.551945760159 
12.519780081038 
21.030654481943 



-10.951590892028 
-10.680835341450 
-9.946495139814 
-8.505293514974 
-6.148266454200 
-2.710246349993 
1.934250402020 
7.879187144780 
15.194895308886 
23.934394836774 
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Table 2. Values of /? as defined in Eq. (|2l|), for Z = 2, 3 < n < 10 , presented as in 
Table I 



n 



4 



6 



-3.278125698868 
-0.798234614319 
4.076360313187 



-4.521965213374 
-2.785523299847 
0.779586310517 
6.527902202704 



-5.668426333247 
-4.411763596407 
-1.735936565934 
2.700152211531 
9.115974284057 



-6.761249358836 
-5.826275313135 
-3.774065245823 
-0.287039774645 
4.853558868056 
11.795070824383 



n = 7 



9 



10 



-7.822716061200 
-7.110457925768 
-5.507083989124 
-2.723146982242 
1.453378478484 
7.172075276697 
14.537951203153 



-8.864931236359 
-8.311162452760 
-7.038000908935 
-4.785217058704 
-1.351453753382 
3.411667214518 
9.611814936391 
17.327283259231 



-9.894830804857 
-9.456572802320 
-8.431024341412 
-6.586194666454 
-3.733658908457 
0.271140940945 
5.536500893854 
12.143135917875 
20.151503770825 



-10.916571653805 
-10.564284465555 
-9.727562574162 
-8.200598472288 
-5.808927517424 
-2.413663893356 
2.092437435636 
7.791260127102 
14.745312392743 
23.002598621108 
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Table 3. Values of f3 as defined in Eq. (pT|), for / = 3, 3 < n < 10, presented as in 
Table I 



n 



4 



6 



-3.151022588503 
-0.688273381129 
3.839295969632 



-4.406309629179 
-2.602022794423 
0.834376786655 
6.173955636947 



-5.569879138144 
-4.214093797042 
-1.552972495811 
2.678798505483 
8.658146925515 



-6.679108747771 
-5.637934713748 
-3.541460362647 
-0.137921230589 
4.750180241001 
11.246244813755 



n = 7 



9 



10 



-7.754622190689 
-6.939787368805 
-5.263146046115 
-2.490737352918 
1.551042361196 
6.987662691921 
13.909587905409 



-8.808385081533 
-8.160018019327 
-6.801150749358 
-4.517404302068 
-1.142156136416 
3.449539586975 
9.350359497993 
16.629215203734 



-9.847641443106 
-9.324116930842 
-8.209581772798 
-6.309419903787 
-3.465931796518 
0.443287284785 
5.511558811618 
11.809781767444 
19.392063983204 



-10.876940137334 
-10.448678615450 
-9.524662708154 
-7.929266399940 
-5.514770174927 
-2.162791993620 
2.219325513867 
7.703479769014 
14.345470691179 
22.188834055364 
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